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Flight test and modeling techniques were developed for efficiently identifying global 
aerodynamic models that can be used to accurately simulate stall, upset, and recovery on 
large transport airplanes. The techniques were developed and validated in a high-fidelity 
fixed-base flight simulator using a wind-tunnel aerodynamic database, realistic sensor 
characteristics, and a realistic flight deck representative of a large transport aircraft. 
Results demonstrated that aerodynamic models for stall, upset, and recovery can be 
identified rapidly and accurately using relatively simple piloted flight test maneuvers. Stall 
maneuver predictions and comparisons of identified aerodynamic models with data from the 
underlying simulation aerodynamic database were used to validate the techniques. 


Nomenclature 


a x ,a y ,a z 

— 

body-axis translational accelerations, g 

h 

= 

wing span, ft 

c 

= 

wing mean aerodynamic chord, ft 

c x ,c Y ,c z 

= 

body-axis nondimensional aerodynamic force coefficients 

C m > C n 

= 

body-axis nondimensional aerodynamic moment coefficients 

I x’ I y’ I z’ I xz 

= 

mass moments of inertia, slug-ft 2 

m 

= 

aircraft mass, slug 

M t 

= 

body-axis pitching moment from engine thrust, ft-lbf 

p, q,r 

= 

body-axis roll, pitch, and yaw rates, rad/s or deg/s 

PSE 

= 

predicted squared error 

q 

= 

dynamic pressure, lbf/ft 2 

s 

= 

wing reference area, ft 2 

S1DPAC 

= 

System IDentification Programs for AirCraft 

T 

= 

maneuver length, s 

X Y , Z Y 

= 

body-axis components of engine thrust, lbf 

V 

= 

true airspeed, ft/s 

a 

= 

angle of attack, rad or deg 

p 

= 

sideslip angle, rad or deg 

S a ,S e ,8 r 

= 

aileron, elevator, and rudder deflections, rad or deg 

(/), 6,\j/ 

= 

Euler roll, pitch, and yaw angles, rad or deg 

I 

= 

covariance matrix 
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superscripts 

T 

-1 

subscripts 

eg 

o 


transpose 
estimate 
time derivative 
matrix inverse 


center of gravity 
reference value or base term 


I. Introduction 

O N August 1, 2010, the Airline Safety and Federal Aviation Administration Extension Act of 2010 (H.R. 5900) 
became law in the United States. This law was intended to improve aviation safety by addressing problems 
that could have played a role in the February 2009 crash of Colgan Air Flight 3407. One of the provisions in the law 
directs the Federal Aviation Administration (FAA) to ensure that commercial airline pilots are trained to recover 
from stalls and upsets, where an upset is generally defined as pitch attitude outside the range [-10, 25] degrees, or 
roll attitude outside the range [-45, 45] degrees, or airspeed inappropriate for the conditions 1,2 . Accurate flight 
simulation for aircraft stall, upset, and recovery is required to enable the effective pilot training necessary to comply 
with the law, and thereby improve aviation safety. 

Typically, flight simulations are based on aerodynamic data collected from wind tunnel tests or generated by 
computational methods such as computational fluid dynamics (CFD). For large transport aircraft, the aerodynamic 
database used for flight simulation often does not include data required for accurate stall upset and recovery 
training 1,2 . Although current flight simulators rely on flight data for validation, flight validation is limited for stall 
upset conditions. The situation has been widely recognized in government and industry, and is being addressed with 
research that expands wind-tunnel testing and CFD to include much larger ranges of motion and control variables 3-8 . 
With expanded flight envelope data, flight simulations can be built to more accurately represent aircraft flight 
characteristics in stall upset and recovery conditions. 

The present work examines an alternate approach for obtaining an accurate aerodynamic model over extended 
flight envelopes. The approach involves efficient piloted flight testing and nonlinear aerodynamic modeling that can 
produce accurate global aerodynamic models based on flight data alone. The advantages of this approach are time 
and cost savings, in addition to better fidelity that can be expected from generating global aerodynamic models 
using direct measurements on the real full-scale airplane in flight, as opposed to using subscale wind-tunnel models 
or mathematical representations of the aircraft and the air flow. The latter approaches have accuracy limitations 
because of effects such as model scale and geometry differences compared to the full-scale aircraft, wind tunnel wall 
and sting interference, wind tunnel flow angularity, Reynolds number differences, turbulence modeling deficiencies, 
grid geometry deficiencies for the aircraft and flow field, and high costs. Using flight test methods to generate the 
aerodynamic model directly avoids all of these problems, although the typical problems associated with any flight 
test still remain, e.g., sensor data quality, executing flight test maneuvers with sufficient data information content, 
high cost, safety risk, and practical considerations. The methods described in this work are intended to be a 
practical, cost-effective, and efficient means to improve flight simulation in stall, upset, and recovery flight 
conditions, to help comply with the law requiring pilot training in these conditions for improved aviation safety. 
Furthermore, the FAA requires flight data matching criteria to be satisfied as part of the certification process for 
commercial flight training simulators. If the aerodynamic model for flight simulation is generated accurately from 
flight data directly, this requirement could be satisfied without further work. 

Previous work in research flight test environments studied effective flight test maneuvers 9 ' 12 and nonlinear 
modeling techniques 11,12 for efficiently generating global aerodynamic models from flight test data alone. The 
present work refines and extends those techniques, and applies them to the practical problem of generating 
high-fidelity flight simulation for large transport airplanes over an expanded flight envelope. The work was done in 
a flight simulation environment, so that the effectiveness of the modeling techniques could be evaluated by directly 
comparing the underlying aerodynamic database used to generate the simulated flight data with the aerodynamic 
models identified from the simulated flight data alone. 

The next section describes the simulation facility used for this work. Following that is a description of the 
piloted flight test maneuver technique, and some evaluation and discussion of the resulting simulated flight data. 
Section IV describes the global aerodynamic modeling technique, which combines multivariate polynomials, 
polynomial splines, orthogonal function modeling theory, and a statistical modeling metric to autonomously identify 
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compact global aerodynamic models for nondimensional force and moment coefficients. The piloted flight test 
maneuvers and global aerodynamic modeling method are then demonstrated using simulated flight data with 
realistic aerodynamic nonlinearities and stall characteristics. Direct comparisons of the underlying aerodynamic 
database and the identified global aerodynamic models are shown in Section V. Prediction results are presented in 
Section VI. 

All of the experiment design, data analysis, and modeling tasks included in this work were done using system 
identification software written in MATLAB®, called System IDentification Programs for AirCraft, or SIDPAC 13 . 
SIDPAC is bundled with Ref. [13], and is therefore publicly available. The SIDPAC software toolbox was 
developed at NASA Langley, and is continually expanded and improved. SIDPAC has been applied successfully to 
a wide variety of flight and wind tunnel experiments at NASA Langley and elsewhere, and is used at more than 80 
institutions worldwide 14 . 


II. Aircraft Simulation 

This study was conducted in the NASA Langley Cockpit Motion Facility using the Integration Flight Deck 
simulator cockpit. Detailed discussions of the facility can be found in Refs. [15] -[16]. The cockpit is shown in 
Fig. 1. The simulation used for this work was a research implementation of a level D crew training flight 
simulation. The airplane simulated was a transcontinental -range commercial jet transport with twin engines 
mounted under the wing and a conventional, low, horizontal tail. Mass/inertia properties were representative of a 
mid-weight condition. For this study, the simulator was operated in a fixed-base mode, which did not use the 
hexapod motion base system. 


A. Flight Dynamics Model and Aerodynamic Database 


B. Flight Controls 

The flight control system used in this simulation 
study is a high-fidelity representation of a 
conventional hydromechanical flight control system 
for a transport airplane. A control column, wheel, and 
rudder pedals were used for primary flight control. 

These control inceptors were connected to a hydraulic 
control loader which accurately simulated the 
characteristics of a force-feel system for a transport 
airplane. The primary control surfaces are elevator, 
stabilizer, left and right ailerons, inboard and outboard left and right spoilers, and rudder. The aerodynamic effects 
of these control surfaces were modeled as increments to the baseline aerodynamic model. For the maneuvers in this 
study, the stabilizer was always kept at the constant initial deflection required for trim. Maneuvers were generally 
conducted with the yaw damper engaged, unless otherwise noted. 


Figure 1. NASA Langley Integration Flight Deck (IFD) 
simulator cockpit 
Credit: NASA Langley Research Center 



The six-degree-of-freedom simulation is a rigid-body model with quasi-static aeroelastic effects and Reynolds 
number scaling effects. Structural modes are not treated as independent states; rather, aeroelastic effects are 
accounted for by a series of dynamic pressure-based corrections that are applied to the rigid-body aerodynamic 
model. Reynolds number effects are incorporated directly into the aerodynamic database. The original crew training 
flight dynamics model was enhanced based on wind 
tunnel measurements to accurately represent the flight 
dynamic characteristic of a transport airplane at stall 
and post-stall flight conditions. The characteristics in 
the aerodynamic database are highly nonlinear with 
respect to relative wind incidence angles and body- 
axis rates at angles of attack near and above stall. 

This enhanced aerodynamic database is described in 
detail in Refs. [3]-[6]. 


C. Simulated Flight Data 

Aircraft sensor data was simulated by applying a sensor measurement noise model to the simulation outputs. 
Systematic instrumentation errors such as scale factor errors and biases were not applied, because typically these 
errors would be removed from measured flight data prior to modeling, using data compatibility analysis 13 . 
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Similarly, simulated aircraft translational accelerations and air flow angles at the aircraft center of gravity were used, 
because flight data would be corrected in this way prior to modeling. Sensor noise levels were specified using 
estimated noise levels obtained by smoothing techniques applied to flight test data from similar maneuvers flown on 
a subscale transport aircraft 10,11,13 . Sensor noise levels are given in Table 1. Sampling rate for the simulated flight 
data was 50 Hz. 


III. Piloted Flight Test Maneuvers 

For local aerodynamic modeling, a pilot would fly the aircraft to a reference flight condition, then execute a local 
perturbation maneuver by applying small perturbation inputs, starting from a steady reference flight condition. For 
global aerodynamic modeling, one conventional approach is to combine local aerodynamic modeling results 
obtained from many of these local perturbation maneuvers at different flight conditions. This requires numerous and 
accurate acquisitions of particular flight conditions, followed by the application of perturbation inputs at each 
condition, and finally combining the local results to produce a global model. 

More efficient global aerodynamic modeling can be achieved by continuously applying multi -axis perturbation 
inputs while the aircraft flight condition is varied slowly. By changing the nominal angle of attack slowly, the 
maneuver can be considered a combination of informative multi-axis excitations for many different nominal angles 
of attack, executed in a single, efficient maneuver. The result is very informative flight data collected efficiently 
over a wide range of nominal angle of attack, which can be used for global aerodynamic modeling. This approach 
was implemented with automated multi-axis excitations while the pilot slowly changed flight condition in subscale 
aircraft flight tests 911 . In the present work, the flight test maneuver was implemented completely by the pilot, using 
an approach analogous to what has been used successfully in the recent past for jet trainer flight tests 12 . 

This section begins with a detailed description of how to fly an efficient flight test maneuver for global 
aerodynamic modeling. Following this, the data from a simulated flight test maneuver flown in this way is 
examined and evaluated for global aerodynamic modeling. 

A. Piloting Technique 

This section describes the piloting technique necessary to execute the flight test maneuvers that generated the 
simulated flight data for global aerodynamic modeling. The general idea is to excite the aircraft dynamics in all 
axes over a short time period by moving multiple control surfaces simultaneously in an uncorrelated manner. The 
excitation inputs should also have wideband frequency content over a range of frequencies encompassing the 
expected modal frequencies for the aircraft dynamic response. These excitation inputs are superimposed on a much 
slower input that slowly changes the flight condition from low angle of attack through stall/upset and recovery. The 
resulting simulated flight data exhibits low correlation among the explanatory variables, over a wide range of flight 
conditions, using a single, efficient piloted maneuver. 

To fly the maneuver, the pilot first established idle thrust trimmed flight at a speed 20 percent over the stall 
speed, which put the airplane in a descent. The stabilizer trim system was configured to maintain a constant 
stabilizer position throughout the maneuver. The pilot then began applying small, uncorrelated perturbation inputs 
to the control column, control wheel, and rudder pedals. At the same time, the pilot slowly applied aft bias to the 
control column as necessary to establish a steady mean pitch attitude of 12.5 degrees, in order to obtain the desired 
airspeed deceleration rate of 2 ft/s per second. Mean pitch attitude was used to maintain a relatively steady airspeed 
deceleration (or steady mean angle of attack increase) because the pitch attitude varied slightly throughout the 
maneuver, as a result of the perturbation inputs. The pilot intentionally varied the frequency of the perturbations and 
made a conscious effort to avoid making correlated inputs, especially between control wheel and rudder pedal. The 
pilot varied the amplitude of the perturbations because stability and control characteristics changed during the 
maneuver as mean angle of attack progressed from low to higher values. The amplitudes of the perturbations were 
sized to obtain perturbation responses in pitch attitude and roll attitude of about ±2 degrees and ±5 degrees, 
respectively. Occasionally, the pilot applied momentary biases to the control wheel perturbations as necessary to 
maintain the desired trajectory and compensate for effects such as stall roll off. 

A good process for implementing the maneuver is to first fly a standard slow approach to stall and recovery, in 
order to establish the feel and airspeed deceleration rate for the slowly -varying flight condition. Then a maneuver 
with similar mean trajectory can be flown, but this time with superimposed uncorrelated perturbation inputs. The 
requirements for a good maneuver are not stringent, so that a wide range of different maneuvers with many 
variations would be considered adequate for the global aerodynamic modeling described later. 
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B. Evaluating the Data from a Piloted Flight Test Maneuver 

Figure 2 shows simulated flight data from a piloted maneuver used to collect data for global aerodynamic 
modeling. Note the simultaneous excitation in all axes, and the slow change in flight condition from low angle of 
attack through high angle of attack, stall, and recovery. 
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Figure 2. Simulated flight data for a global aerodynamic modeling maneuver 

Figure 3 shows cross-plots of aircraft states and controls using data from the maneuver shown in Fig. 2. These 
plots demonstrate that a wide range of the explanatory variables generally used for aerodynamic modeling was 
swept through during this single maneuver. High (undesirable) pair-wise correlations would be indicated by lines or 
narrow ellipses in these plots, whereas nondescript shapes or wide ellipses, as seen in Fig. 3, indicate low (desirable) 
pair-wise correlations. Low pair-wise correlations mean that the aerodynamic dependencies on the explanatory 
variables can be identified accurately and without ambiguity. Cross -plots for other aircraft states and controls were 
similar in that the plots indicated low pair-wise correlations for the explanatory variables. 

This novel maneuver exhibits a combination of low pair-wise correlations for the explanatory variables, multi- 
axis excitation, and slowly-varying flight condition covering large ranges of the explanatory variables to be used for 
aerodynamic modeling. These characteristics make the maneuver very effective and efficient for global 
aerodynamic modeling. 

A series of these piloted maneuvers were flown on the fixed-base commercial jet transport aircraft simulator 
described in Section II. Each maneuver had a constant idle power setting and fixed aircraft configuration (flaps up, 
landing gear up), and all used multi-axis perturbation inputs implemented by the pilot, similar to what is shown in 
Figs. 2 and 3. 
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Figure 3. Simulated flight data cross plots for a global aerodynamic modeling maneuver 


IV. Global Aerodynamic Modeling 

The objective of global aerodynamic modeling is to identify a model for each nondimensional aerodynamic force 
and moment coefficient as a function of explanatory variables that can be measured, such as angle of attack, pitch 
rate, and control surface deflections, over a large range of the explanatory variables. This requires two main tasks: 
1) experiment design to collect dynamic data over a large range of the explanatory variables, which was addressed in 
Section III, and 2) identifying an accurate global model, which is the subject of this section. 

In the typical approach to global aerodynamic modeling, a series of linear or simplified local models, which are 
valid for relatively small ranges of the explanatory variables, are joined together in some way to implement a global 
aerodynamic model. In another approach, a multivariate orthogonal function modeling technique can be used to 
identify global models, which are generally nonlinear and valid over relatively large ranges of the explanatory 
variables 17 . However, for very large ranges of the explanatory variables, or severe local nonlinearities, this approach 
in some cases results in compromised local model fit in order to achieve a better global model fit. This is the result 
of using global orthogonal polynomial functions for the modeling, which do not have the ability to change locally 
without modifying the entire model. Other approaches have used data partitioning with simplified local model 
structures 18 and localized modeling functions 19 to address this problem, and also to make local model updates easier. 
However, these approaches require a procedure or analyst judgment for: 1) partitioning the explanatory variable 
space; 2) selecting the location of the local modeling functions in the explanatory variable space; and/or 3) selecting 
the mathematical structure of the local models. 

In this work, the multivariate orthogonal function modeling developed in Ref. [17] was applied using both the 
explanatory variables and spline functions of the explanatory variables, as described in Ref. [11]. This approach 
provides the required local nonlinear modeling capability, while retaining easy physical interpretation of the model, 
and automating the selection of optimal modeling complexity necessary to accurately characterize the functional 
dependencies. 

The next sections describe the approach to global aerodynamic modeling using multivariate orthogonal functions 
derived from multivariate polynomials and spline functions computed using measured explanatory variable data. 
Further details on this approach to global aerodynamic modeling can be found in Ref. [11]. 

A. Aircraft Aerodynamic Modeling 

For global aerodynamic modeling from flight data, nondimensional aerodynamic force and moment coefficients 
are used as the response variable (also called the dependent variable) in the modeling problem. A separate modeling 
problem is solved for each force or moment coefficient, corresponding to minimizing the equation error in each 
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individual equation of motion for the six rigid-body degrees of freedom of the aircraft. Values for the 
nondimensional aerodynamic force and moment coefficients cannot be measured directly in flight, but instead must 
be computed from measured and known quantities using the following equations 13 



These equations retain the full nonlinear dynamics in the aircraft equations of motion. 

The result is N values of the nondimensional force and moment coefficients, where N is the number of data 
points. These values are often called measured force and moment coefficients, even though they not measured 
directly, but rather computed from other measurements and known quantities. Mass/inertia properties are assumed 
to be known as a function of fuel state. Engine thrust and geometric properties such as wing span and wing 
reference area are assumed to be known. Dynamic pressure, translational accelerations, and explanatory variables 
such as angle of attack, sideslip angle, pitch rate, and control surface deflections, are measured directly. Angular 
accelerations are obtained by smoothed numerical differentiation of the measured angular rates 13 . 

The desired form of the global aerodynamic model is a mathematical model structure with estimated model 
parameter values and associated uncertainties, relating the nondimensional aerodynamic force and moment 
coefficients to aircraft states and controls that can be measured. 

B. Multivariate Orthogonal Function Modeling 

The form of a multivariate orthogonal function model is 

z=a lPl +a 2 p 2 + ... + a nPn +e (6) 

where z is an A-dimensional vector of the response variable (e.g., nondimensional force or moment coefficient), 
z = [ Zi ,z 2 ,-,Zn] T , modeled in terms of a linear combination of n mutually orthogonal modeling functions 
Pj,j = l,2,...,n. Each pj is an A-dimensional vector which in general depends on the explanatory variables. The 
cij , j = 1,2,..., n are constant model parameters to be determined, and £ denotes the modeling error vector. 

Equation (6) represents a mathematical model used to represent functional dependencies in the measured data. 
The important questions of determining how the modeling functions p ; are computed from the explanatory 

variables, as well as which modeling functions should be included in Eq. (6), which implicitly determines n , will be 
addressed later. At this point, the properties of a multivariate orthogonal function model are examined. 

Define an Nxn matrix P , 

P = [PvP2>-’Pn] ( 7 ) 
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and let a - [a 1 ,a 2 ,.~,a n T . Equation (6) can then be written as a standard least squares regression problem, 


z = Pa + s (8) 

The error vector s is to be minimized in a least squares sense. The goal is to determine a that minimizes the least 
squares cost function 

J=-^£ T £ = -^(z-Pa) T (z- Pa) (9) 

The parameter vector estimate a that minimizes this cost function is 13 

a = [ pT p]~ lpT z ( 10 ) 

The estimated parameter covariance matrix is 13 




(d-a)(d-af = a 2 (p T p) _1 


(ID 


where E is the expectation operator, and the fit error variance <j z can be estimated from the residuals 

v =z~ Pci (12) 


using 


ci 


1 

(N-n) 


( z-Paf ( z-Pa ) 


(N-n) 


(13) 


Parameter standard errors are computed as the square root of the diagonal elements of the 2^ matrix from Eq. (11), 

^ 9 

using cr from Eq. (13). The identified model output y is 


y = Pa 


(14) 


In conventional least-squares modeling, the modeling functions (columns of P) are often polynomials in the 
explanatory variables. This approach corresponds to using the terms of a multivariate Taylor series expansion to 
approximate the functional dependence of the response variable on the explanatory variables. If the modeling 
functions are instead multivariate orthogonal functions generated from the explanatory variable data, it is easier to 
determine an appropriate model structure, because the explanatory capability of each modeling function is 
completely distinct from all the others. This decouples the least squares modeling problem, as will be shown now. 

For mutually orthogonal modeling functions, 


Pi Pj =0 , i*j , i, j = 1, 2 , n 


(15) 


and P T P is a diagonal matrix with the inner product of the orthogonal functions on the main diagonal. Using 
Eqs. (7) and (15) in Eq. (10), the / h element of the estimated parameter vector a is computed as 

a j=(pjz)/(pj Pj) < 16 > 
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Using Eqs. (7), (15), and (16) in Eq. (9), 





( 17 ) 


Equation (17) shows that when the modeling functions are orthogonal, the reduction in the least-squares cost 
function resulting from including the term aj pj in the model depends only on the response variable data z and the 

added orthogonal modeling function pj . The least-squares modeling problem is therefore decoupled, which means 

each orthogonal modeling function can be evaluated independently in terms of its ability to reduce the least- squares 
model fit to the data, regardless of which other orthogonal modeling functions are already selected for the model. 
When the modeling functions p j are instead polynomials in the explanatory variables (or any other non -orthogonal 

function set), the least-squares problem is coupled, and iterative analysis is required to select the modeling functions 
for an adequate model structure. 

The orthogonal modeling functions to be included in the model are chosen to minimize predicted squared error, 
PSE , defined by 20 


re£ . ( Z -P«) (z-K) + j n_ 

N max n 


(18) 


or 


dce ,_ 2 / 2 n 

PSE — h cr max — 

N N 


(19) 


The constant cj} nax is the upper-bound estimate of the squared error between future data and the model, i.e., the 
upper-bound mean squared error for prediction cases. The upper bound is used in the model over-fit penalty term to 
account for the fact that PSE is calculated during the model structure determination stage, when the model structure 

is not correct. Using <J^ nax in the model complexity term also results in the PSE metric conservatively estimating 
the squared error for prediction cases. 

A simple estimate of cr^ ax that is independent of the model structure can be obtained by considering <t^ x to be 
the residual variance estimate for a constant model equal to the mean of the measured response values, 


N - 1 


N 

Zh- ■ 

i=\ 


where 



( 20 ) 


( 21 ) 
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The PSE in Eq. (19) depends on the mean squared error, 2j/N , and a term proportional to the number of terms 
in the model, n . The latter term prevents over-fitting the data with too many model terms, which is detrimental to 
model prediction accuracy 13,20,21 . While the mean squared error 2 j/N must decrease with the addition of each 
orthogonal modeling function to the model (by Eq. 

(17)), the over-fit penalty term <J^ nax n/N must increase 
with each added model term (increased n). Introducing 
the orthogonal modeling functions into the model in 
order of most effective to least effective in reducing the 
mean squared error (quantified by 

(2/ TV ) [p T jZ j j^p 1 - Pj^ for the / h orthogonal modeling 

function) results in the PSE metric always having a 
single global minimum. 

Figure 4 depicts this graphically, using actual 
modeling results from Ref. [21]. The figure shows that 
after the first 6 modeling functions, the added model 
complexity associated with an additional orthogonal 
modeling function is not justified by the associated 
reduction in mean squared error. This point is marked 
by minimum PSE, which defines an adequate model 
structure with good predictive capability. Ref. [20] 
contains further statistical arguments and analysis for the form of PSE given in Eq. (18), including justification for 
its use in modeling problems. 

Using orthogonal functions to model the response variable makes it possible to evaluate the merit of including 
each modeling function individually , using the predicted squared error PSE. The goal is to select a model structure 
with minimum PSE, and the PSE always has a single global minimum for orthogonal modeling functions. This 
makes the model structure determination a well-defined and straightforward process that can be (and was) 
automated. 

C. Generating Orthogonal Modeling Functions 

Multivariate orthogonal functions can be generated from ordinary multivariate functions in the explanatory 
variables using a Gram-Schmidt orthogonalization procedure. This approach is described in Refs. [11], [13], [17], 
and [21], which are the basis for the material presented here. 

The process begins by choosing one of the ordinary multivariate functions as the first orthogonal function. 
Typically, a vector of ones (associated with the bias term in the model) is chosen as the first orthogonal function, 

Pi — 1 (22) 

In general, any function of the explanatory variables can be chosen as the first orthogonal function, without any 
change in the procedure. To generate the next orthogonal function, an ordinary multivariate polynomial function is 
made orthogonal to the preceding orthogonal function(s). Define the / h orthogonal function pj as 


x 10 5 



Orthogonal Function Number 


Figure 4. Model structure determination using 
orthogonal functions and PSE 


j - 1 

Pj = £/ - X fly Pk J = 2, 3, n, (23) 

k= 1 

where £ j is the / h ordinary multivariate function vector. For example, each could be some ordinary polynomial 

function of the explanatory variables, or spline functions of the explanatory variables, or multiplicative 
combinations of these, or any general multivariate function. The y k j for k = 1, 2, — 1 are scalars determined by 

multiplying both sides of Eq. (23) by p\ , then invoking the mutual orthogonality of the p k , k = and 

solving for y k j , 
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* = 1,2,..., ;-l 


(24) 


ru = 


Pktj 

PkPk 


The same process can be implemented in sequence for each ordinary multivariate function €j,j = 2,3,...,n t . 

The total number of ordinary multivariate functions used as raw material for generating the multivariate orthogonal 
functions, including the bias term, is n t . It can be seen from Eqs. (22)-(24) that each orthogonal function can be 
expressed exactly in terms of a linear expansion of the original multivariate functions. The orthogonal functions are 
generated sequentially by orthogonalizing the original multivariate functions with respect to the orthogonal 
functions already computed. Each orthogonal function can be considered an orthogonalized version of an original 
multivariate function. 

The multivariate orthogonal function generation method described here normally starts by generating all possible 
ordinary multivariate polynomials in the explanatory variables, up to a selected maximum order. For example, when 
modeling the nondimensional vertical force coefficient C z , the explanatory variables might be angle of attack a , 
nondimensional pitch rate q = qc/ 2V, and elevator deflection S e . If the selected maximum order is 3, then the 
ordinary multivariate polynomial modeling functions used as the raw material for the orthogonalization process 
would include terms like a,q,S e ,a ,S e ,aq ,aqS e , etc. Note that considering any other candidate explanatory 
variables, such as sideslip angle p , can be done by simply including the sideslip angle among the group of 
explanatory variables. If it turns out that the sideslip angle is not needed to model C z , the model structure 
determination using orthogonal modeling functions will not select any orthogonal functions associated with sideslip 
angle. This occurs automatically in the course of the model structure determination process described earlier. 
Therefore, there is no harm in including explanatory variables that might not be important, except that additional 
computation time will be required to identify the model structure, because additional multivariate orthogonal 
functions will be generated and sorted. Similarly, if the maximum order is chosen higher than necessary, the only 
penalty would be the increased computation time necessary for generating and sorting the additional orthogonal 
functions. The final identified model would be the same. Consequently, the choices that the analyst needs to make 
are easy and not critical to the quality of the final modeling results. 

If the pj vectors and the gj vectors are arranged as columns of matrices P and X , respectively, and the y k j 
are elements in the k th row and / h column of an upper triangular matrix G with ones on the diagonal, 



"1 

Yn 

y 13 ■ 

y 1 n t 


0 

1 

y 23 ■ 

" Yin, 

G - 

0 

0 

i • 

■■ Y3n, 


0 

0 

o • 

•• 1 


Then 

X = PG 


which leads to 


P = XG~ l 


(25) 


(26) 


(27) 


The columns of G 1 contain the coefficients for expansion of each column of P (i.e., each multivariate orthogonal 
function) in terms of an exact linear expansion in the original multivariate functions in the columns of X . 
Equation (27) can be used to express each multivariate orthogonal function in terms of the original multivariate 
functions. The manner in which the orthogonal functions are generated allows them to be decomposed without 
ambiguity into an expansion of the original multivariate functions, which have physical meaning. 
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D. Conversion to Physically Meaningful Multivariate Function Models 

After the model structure is determined using multivariate orthogonal modeling functions for minimum PSE , the 
identified model output is 


y=Pd (28) 

where the P matrix now includes only the n orthogonal functions selected in the model structure determination, 
n<n t . Each retained orthogonal modeling function can be decomposed without error into an expansion of the 

original multivariate functions in the explanatory variables, using the columns of G~ l in Eq. (27) corresponding to 
the retained orthogonal functions. Common terms are combined using double precision arithmetic to arrive finally 
at a model using only original multivariate functions in the explanatory variables. Terms that contribute less than 
0.1 percent of the final model root-mean- square magnitude are dropped. 

The final form of the model is a sum of ordinary multivariate functions in the explanatory variables, with 
associated model parameter estimates and uncertainties. 

E. Including Spline Functions in the Orthogonalization 

For global aerodynamic modeling, functional dependencies of the aerodynamic coefficients on the explanatory 
variables can exhibit significant localized variation or changes in character. In these cases, a global polynomial 
model can be inadequate for capturing the local variations. To solve this problem, spline functions in the 
explanatory variables can be introduced as additional explanatory variables. Splines have local modeling capability 
that global polynomial modeling functions do not have, because splines are polynomials defined only on selected 
intervals. Low-order polynomial terms defined on limited intervals can approximate local nonlinearities quite well. 
Using polynomial splines also retains clear physical interpretation in the final model. 

Spline functions are defined as piecewise polynomials functions of degree m in one or more explanatory 
variables. The term “piecewise” means that the polynomial is different for specific ranges of the explanatory 
variables. Spline function values and derivatives agree at the points where the piecewise polynomials join. These 
points are called knots, and are defined as specific values of each explanatory variable. A polynomial spline S m (x) 

of degree m with continuous derivatives up to degree m-1 , for a single explanatory variable x e \xrmn> x max\ > can 
be expressed as 


where 


S m (*) = Z C r X " + E D i ( X ~ X i )- 

r = 1 i = 1 


/ \m \(x-X i ) n 

(*-*)+ =j V ” 


X > x t 
x<x t 


(29) 


(30) 


and the C r and D t are constants and are the piecewise parts of the spline function. The values 

x lf x 2 , . . x k are knots which satisfy the condition 

x rmn< x \< x 2<---< x k< x mwc ( 31 ) 

Three piecewise spline types are sketched in Fig. 5 for the same selected knots in angle of attack. 

Note that when the spline knots are the same, a higher-order piecewise spline in a single explanatory variable can 
be computed as a multiplication of lower-order piecewise splines in the same explanatory variable, 

={ x - x i)+~ 1 ( x ~ x if+ m - 1 ( 32 ) 
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This fact can be used advantageously in generating multivariate orthogonal modeling functions with excellent 
local modeling capability. An example will be used to demonstrate the idea. 

Suppose that the nondimensional vertical aerodynamic force coefficient C z is being modeled with explanatory 
variables angle of attack a and elevator control deflection 8 e . For a maximum selected model order 2, the ordinary 
polynomial modeling functions that will serve as raw material for the orthogonalization process are: 

1, a , a 2 , aS e , 8 e , S 2 (33) 


which would lead to a multivariate polynomial model of the form 


c z ~ c z n + c z„ a + c z 9 Ct 2 + C z aS e +C z S e +c z S 2 


(34) 


where the values of the model parameters such as C z ^ and C Zag would be estimated from the data using the 

orthogonal function modeling and subsequent decomposition procedures described earlier. Some of the terms on the 
right side of Eq. (34) might not be present, depending on the results of the model structure determination using 
orthogonal functions. This represents a global polynomial model. 

Now introduce a first-order spline in angle of attack with a single knot located at 10 deg as an additional 
explanatory variable. Using the same maximum model order 2, the set of ordinary polynomial modeling functions 
expands to 

1, a, a 2 , a8 e , 8 e , S 2 ,(a-10)\ , a(a-l0) l + , 8 e (er-10)^ , (a-l0) 2 + (35) 


which would lead to a multivariate polynomial model of the form 

c z= c z„+ c z a a + c z A + Cz as e + C z S e + C z S 2 


+ C Zi (a-lO)\+C z r a(a-lO)[ + C z f S e (a-10)[ +C Z 2 (or- 10)* 


(36) 


As before, some of the terms on the right side of Eq. (36) 
might be not be present, depending on the results of the 
model structure determination using orthogonal functions. 

With this expanded set of ordinary modeling 
functions as raw material, more orthogonal functions will 
be generated, and there is now capability to accommodate 
changes in the linear dependence of C z on a for 

a >10 deg , as well as additional a and aS e 

nonlinearities that only take effect when a > 10 deg . 

Note also that a(a- 10)^ and (<z-10)^, for example, 

are in general not the same function. Because 
higher-order splines are created automatically by the 
multiplication of first-order splines (cf. Eq. (32)), only 
first-order splines in the explanatory variables need to be 
included as additional explanatory variables. 

Extrapolating from this simple example, it is clear that 
multiple knots in the explanatory variables would provide 
local nonlinear modeling capability in multiple 
explanatory variables, while retaining physical insight 
into the functional dependencies. The automated 
orthogonalization and sorting process described earlier 
identifies which nonlinear terms are necessary to 



a (deg) 

Figure 5. Polynomial splines in angle of attack: 
(a) zero-degree, (b) first-degree, (c) second-degree 
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characterize the functional dependencies and estimates the associated model parameters and associated uncertainties. 
Inputs required from the analyst relate only to the limits of what should be considered, such as which explanatory 
variables to consider, maximum model order, and knot locations. However, these can be specified very generously, 
because the orthogonal function modeling algorithm automatically sorts out which of the terms are important, based 
on the data, and discards the rest. The result is a global parsimonious model with excellent global and local 
nonlinear modeling capability and easy physical interpretation. 

When many different splines are considered for closely-spaced knot locations, perhaps for more than one 
explanatory variable, it is possible to run into memory limitations for the large number of orthogonal functions that 
are generated. This problem can be avoided by considering subsets of the full set of splines during a series of 
identification stages, then retaining only the selected spline terms at successive stages. This amounts to a sequential 
search for effective modeling terms. For example, the first stage might consider only first-order and second-order 
splines in angle of attack, with knots spaced 1 deg apart, from 8 deg to 22 deg. The model structure determination 
might then select only first-order spline terms at angle of attack knot locations of 16 deg and 20 deg. The next stage 
would retain all of the terms identified for the model in the first stage, but might investigate first-order and second- 
order splines in sideslip angle from -10 deg to 10 deg with knots spaced 1 deg apart, for example, and so on. The 
general idea is that once a candidate nonlinear or spline modeling function is discarded at any stage, it need not be 
considered again. This saves computation and memory during the model structure identification process, while still 
allowing a wide range of nonlinear and spline modeling terms to be considered for inclusion in the final model. 

V. Results 

The techniques for piloted flight test maneuvers and global aerodynamic modeling described in Sections III and 
IV were applied using the fixed-base commercial jet transport aircraft simulator described in Section II. Data from a 
typical piloted maneuver is shown in Figs. 2 and 3. 

A. Identified Models and Comparison with the Simulation Database 

Applying the flight testing and global modeling techniques in simulation allows comparison of the global 
aerodynamic models identified from the simulated flight test data with the known underlying aerodynamic database 
used in the simulation. The aerodynamic database in the simulation was implemented in tabular form, based on 
wind tunnel data. Simulated flight data from the piloted maneuver shown in Figures 2 and 3 was used initially for 
the global aerodynamic modeling. Results showed accurate characterization at low angles of attack, but inadequate 
characterization of the relatively complex functional dependencies near stall. It was found that the modeling near 
stall could be improved if an additional stall maneuver without pilot excitation was included in the data used for 
global modeling. The additional maneuver provided enhanced data information content for a more detailed and 
accurate global model near stall. Figure 6 shows the stall maneuver that was combined with the data from Figs. 2 
and 3 to generate the modeling results shown in this section. 

For the global aerodynamic modeling, initial modeling was done using coarse spacing of knot locations in each 
explanatory variable to identify which explanatory variables required the increased local nonlinear modeling 
capability provided by spline terms. This analysis showed that the global aerodynamic modeling required spline 
terms mainly for angle of attack. First-order splines in angle of attack with knots finely spaced at 1 deg intervals 
between 12 deg and 22 deg inclusive (knots = [12:1:22] deg) were added as additional explanatory variables for the 
global aerodynamic modeling. These terms were culled from a larger candidate pool of first-order splines in angle 
of attack with more knot locations, using the automated orthogonal function modeling process. For each 
nondimensional force and moment coefficient model, model terms with up to second-order complexity were 
orthogonalized and sorted to identify an adequate model structure and estimate the associated model parameter 
values and uncertainties. This upper-bound second-order complexity applied to the conventional explanatory 
variables (such as angle of attack, nondimensional pitch rate, and elevator control deflection), as well as the first- 
order splines in angle of attack, which were included as additional explanatory variables. The second -order 
complexity automatically included up to second-order multiplicative combinations of these functions. Because the 
first-order spline terms were multiplied by themselves and each other, up to second order, this approach allowed for 
up to second-order spline terms in angle of attack. Other modeling terms considered in the model identification 

process included nonlinear cross terms invoked at specific values of the angle of attack knots, such as (a-10) 1 ^ S e 

or (or- 12)+ [qc / 2V) . As expected, only some of all the candidate model terms were selected for inclusion in the 
global model for each aerodynamic coefficient. 
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Figure 6. Simulated flight data for a stall maneuver 

Subsequent to the angle of attack spline term identification, a first-order spline in sideslip angle with knot 
location at 0 deg, and zero-order splines in angle of attack with knot locations at 12.4, 12.6, and 12.8 deg were 
identified as important model terms for some of the aerodynamic coefficient models. These additional spline terms 
were identified in a sequential fashion, as described earlier, and were selected from a pool of candidates. It was 
found that additional candidate model terms could be supplied liberally and sequentially to the automated orthogonal 
function modeling procedure, with the result that only important and statistically significant model terms would be 
selected and retained in the model. 

The overall approach permits simple specification of allowable model complexity that results in a remarkably 
diverse and powerful set of candidate modeling terms, from which an adequate model structure is selected 
automatically based on data information content, orthogonalization, and statistical modeling metrics, as explained in 
Section IV. 

Figure 7 shows three-dimensional plots for all six aerodynamic coefficients, comparing the simulated flight data 
(x markers) with a mesh surface drawn by interrogating the aerodynamic database in the simulation, and a smooth 
surface drawn using the global aerodynamic model identified from simulated flight data alone. The ranges of the 
explanatory variables used to generate the mesh surface and the smooth surface were chosen to include all of the 
simulated flight data. Note that numbers on the vertical scales had to be removed because of proprietary restrictions 
on the wind-tunnel database. The results in Fig. 7 show that for all aerodynamic coefficients, the smooth surface 
and mesh surface are very similar, indicating that the global aerodynamic modeling approach captured the 
underlying functional dependencies accurately and efficiently, for all six aerodynamic coefficients simultaneously, 
using data from just two flight test maneuvers. The modeling technique successfully identified complex local 
nonlinearity, as shown in Fig. 7. 

Close inspection shows only small mismatches between the smooth surface and the mesh surface at locations 
relatively far from the simulated flight data, which would be areas of extrapolation for the identified global model. 
Flying another maneuver through those parts of the explanatory variable subspace would provide the data necessary 
to correct the small mismatches. The three-dimensional plots shown in Fig. 7 must necessarily omit some of the 
explanatory variables (because there are more than two explanatory variables, and the third dimension is used for the 
response variable), so the mismatch between simulated flight data and the smooth and mesh surfaces represents 
mostly dependence on explanatory variables that cannot be shown on a three-dimensional plot. 
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Figure 7. Global aerodynamic models identified from simulated flight data 

The results shown in Fig. 7 suggest that a few well-flown maneuvers of the type described here could very 
effectively cover a large portion of the explanatory variable subspace, and when combined with the orthogonal 
function modeling procedure, could be used in place of wind-tunnel testing that requires extensive grids of test 
points to generate the aerodynamic database. Because the modeling uses an equation -error least- squares 
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formulation, identifying a model based on combined data from multiple maneuvers is straightforward, and is done 
by simply stacking the data from multiple maneuvers into a single larger problem 22 . 

As an example of the compact and analytical final form of the global model, the global model identified for 
nondimensional aerodynamic vertical force coefficient C z from simulated flight data was 

C z =C Zo +C Z( a + C Zi (a-13f ++ C Zo (a-12.4)° +C z<a aS e +C Z , (cc-U)\ 


-C Zo (a-12.Z)l+C Z] ( a -20)\ + C z a^- + C Zi (a-l6)\ + C Zg S e 


(37) 


+ C ^f( B - 12 i +C ^f + S (B - 19t 


where model parameter estimates and uncertainties cannot be given, because of proprietary data restrictions. 
Although data restrictions significantly compromised the presentation of results in this work, the fact that the wind- 
tunnel database is considered highly realistic 3 ' 6 was judged to be more important. Showing that the piloted 
maneuvers and global modeling techniques can be successfully applied in simulation using this realistic database is 
a strong indicator of potential success for application to real flight testing of large commercial jet transport aircraft. 

B. Prediction Testing 

The prediction capability of the identified global aerodynamic models was tested using simulated flight data 
from maneuvers that were not used in the global model identification. Figure 8 shows data from a simulated slow 
approach to stall without pilot excitations, with the yaw damper on. The predictions shown in Fig. 8 are remarkably 
good, considering that data from only two flight test maneuvers was used to identify the global models. Figure 9 
shows the prediction for a similar slow approach to stall maneuver, but with the yaw damper turned off. In this case, 
there is some degradation in the prediction, particularly for the yawing moment coefficient near stall. This is likely 
because the explanatory variables for this maneuver (where yaw damper was off) ranged beyond the 
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Figure 8. Stall maneuver prediction, yaw damper on 
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subspace covered by the explanatory variables for the identification data (where yaw damper was on), making this 
prediction case at least a partial extrapolation for the lateral coefficients. Note that the quality of the longitudinal 
coefficient predictions in Fig. 9 was similar to that shown in Fig. 8, because the yaw damping affects mainly lateral 
explanatory variables. For both of the maneuvers shown in Figs. 8 and 9, the root-mean-square prediction error was 
less than 0.01 for the longitudinal coefficients (C x ,C z ,C m ) , and less than 0.001 for the lateral coefficients 

(Cy , Q, C n ) . Note that these prediction errors included the measurement noise. 
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Figure 9 . Stall maneuver prediction, yaw damper off 

VI. Conclusions 

Simulated flight test data for a commercial jet transport aircraft was used to demonstrate techniques for 
efficiently and effectively identifying global aerodynamic models from piloted flight test maneuvers. The identified 
global aerodynamic models were shown to be accurate, with good prediction capability and easy physical 
interpretation. The identified models can be used for accurate flight simulation of stall, upset, and recovery, as well 
as normal operational flight. 

A piloted flight test maneuver technique was developed and used to collect data over relatively large ranges of 
the explanatory variables for all six rigid-body degrees of freedom of the aircraft simultaneously. This resulted in 
efficient collection of highly informative dynamic flight data with de -correlated explanatory variables, for good 
model identifiability and modeling accuracy. The global aerodynamic modeling method incorporated spline 
functions in an automated orthogonal function modeling scheme, to allow accurate modeling of local complicated 
functional dependencies without adverse effects on other parts of the global model. 

The efficiency of the approach came from the fact that the global aerodynamic modeling was automated, based 
on flight data alone, and also from the efficiency of the piloted flight test maneuvers. The required maneuvering is 
not difficult to learn or implement, and does not require flying the aircraft in unusual or dangerous flight conditions 
for long periods of time. The resulting global aerodynamic models exhibited good prediction capability for data 
from flight maneuvers that were not used in the modeling process. Global aerodynamic models identified from 
flight data could also serve as the basis for certifying flight simulator fidelity, e.g., using Federal Aviation 
Administration Level D certification standards for flight simulators. 
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In addition to the intended use for accurate and efficient stall, upset, and recovery modeling, the experimentation 
and modeling techniques developed in this work could be used in rapid and relatively low-budget unmanned aircraft 
programs to save development time and money by making it possible to generate an accurate nonlinear aircraft 
simulation from a single flight, without the need for extensive wind tunnel testing or aerodynamic calculations. A 
global aerodynamic model identified from flight data alone could also be compared to aerodynamic databases 
generated from wind tunnel data and aerodynamic calculations to improve the fidelity of those ground-based 
aerodynamic prediction methods. This could be done by identifying where there are significant differences between 
the database generated using ground-based methods and the global model identified from flight data alone, as well 
as where the models agreed well. Such comparisons could also be used to estimate full-scale Reynolds number 
effects, or artifacts of wind tunnel testing, such as sting interference and wall effects. 

The techniques described here could be applied to other commercial aircraft, such as regional jets and turboprop 
aircraft, and also to general aviation or military aircraft. The improved efficiency associated with this flight testing 
approach makes the idea of identifying an accurate flight simulation from flight data for each individual aircraft (as 
opposed to one simulation for all aircraft of the same make and model) a reasonable proposition. This capability 
could be a useful tool for improving aircraft maintenance, fault detection, envelope protection, and aviation safety. 

VII. Acknowledgments 

This research in Aircraft System Identification was funded by the NASA Aviation Safety Program, Vehicle 
Systems Safety Technologies (VS ST) project. 

VIII. References 

^ambregts, A.A., Nesemeier, G., Wilbom, J.E., and Newman, R.L., “Airplane Upsets: Old Problem, New Issues”, AIAA 
Modeling and Simulation Technologies Conference and Exhibit , AIAA 2008-6867, Honolulu, HI, August 2008. 

2 Advani, S.K., Schroeder, J.A., and Burks, B., “What Really Can Be Done in Simulation to Improve Upset Training AIAA 
2010-7791, AIAA Modeling and Simulation Technologies Conference, Toronto, Ontario Canada, August 2010. 

3 Foster, J.V., Cunningham, K., Fremaux, C.M., Shah, G.H., Stewart, E.C., Rivers, R.A., Wilborn, J.E., and Gato, W. 
“Dynamics Modeling and Simulation of Earge Transport Airplanes in Upset Conditions,” AIAA 2005-5933, AIAA Guidance , 
Navigation, and Control Conference and Exhibit, San Francisco, CA, August 2005. 

4 Shah, G.H., Cunningham, K., Foster, J.V., Fremaux, C.M., Stewart, E.C., Wilbom, J.E., Gato, W., Pratt, D.W., “Wind- 
Tunnel Investigation of Commercial Transport Aircraft Aerodynamics at Extreme Flight Conditions”, SAE 2002-01-2912, World 
Aviation Congress & Display, Phoenix, AZ, November 2002. 

5 Cunningham, K., Foster, J.V., Shah, G.H., Stewart, E.C., Rivers, R.A., Wilbom, J.E., and Gato, W., “Simulation Study of a 
Commercial Transport Airplane During Stall and Post-Stall Flight”, SAE 2004-01-3100, 2004 SAE World Aviation Congress & 
Exposition, Reno, NV, November 2004. 

6 Cunningham, K., Foster, J.V., Shah, G.H., Stewart, E.C., Ventura, R.N., Rivers, R.A., Wilborn, J.E., and Gato, W., 
“Simulation Study of Flap Effects on a Commercial Transport Airplane in Upset Conditions”, AIAA 2005-5908, AIAA 
Atmospheric Flight Mechanics Conference and Exhibit, San Francisco, CA, August 2005. 

7 Murch, A.M. and Foster, J.V. “Recent NASA Research on Aerodynamic Modeling of Post-Stall and Spin Dynamics of 
Earge Transport Airplanes”, AIAA-2007-463-253, 45 th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, January 2007. 

8 Shah, G.H., Foster, J.V., Cunningham, K., and Owens, D.B. “Simulation Modeling for Off-Nominal Conditions - Where 
Are We Today?”, AIAA 2010-7792, AIAA Modeling and Simulation Technologies Conference, Toronto, Ontario Canada, August 
2010 . 

9 Cunningham, K., Foster, J.V., Morelli, E.A., and Murch, A.M., “Practical Application of a Subscale Transport Aircraft for 
Flight Research in Control Upset and Failure Conditions,” AIAA-2008-6200, AIAA Atmospheric Flight Mechanics Conference, 
Honolulu, HI, August 2008. 

10 Morelli, E.A. “Flight Test Maneuver Design for Efficient Aerodynamic Modeling,” AIAA 2011-6672, AIAA Atmospheric 
Flight Mechanics Conference, Portland, OR, August 2011. 

1 Corelli, E.A. “Efficient Global Aerodynamic Modeling from Flight Data,” AIAA-2012-1050, 50 th AIAA Aerospace 
Sciences Meeting, Nashville, TN, January 2012. 

12 Brandon, J.M. and Morelli, E.A. “Nonlinear Aerodynamic Modeling From Flight Data Using Advanced Piloted Maneuvers 
and Fuzzy Fogic ,” NASA/TM-2012-217778, October 2012. 

13 Klein, V. and Morelli, E.A., Aircraft System Identification - Theory and Practice, AIAA Education Series, AIAA, 
Reston, VA, 2006. 

14 http://dcb. larc.nasa.gov/SIDPAC/SIDBook_SIDPAC.htm 

15 Smith, R.M.: “A Description of the Cockpit Motion Facility and the Research Flight Deck Simulator,” AIAA-2000-4174, 
AIAA Modeling and Simulation Technologies Conference and Exhibit, Denver, CO, August 2000. 

16 Chung, V. and Hutchinson, B.: “A Unique Software System for Simulation-To-Flight Research,” AIAA-200 1-4057, AIAA 
Modeling and Simulation Technologies Conference and Exhibit, Montreal, Quebec, August 2001. 


19 

American Institute of Aeronautics and Astronautics 



17 Morelli, E.A. “Global Nonlinear Aerodynamic Modeling using Multivariate Orthogonal Functions,” Journal of Aircraft, 
Vol. 32, No. 2, March-April 1995, pp. 270-77. 

18 Seher-Weiss, S. “Identification of Nonlinear Aerodynamic Derivatives using Classical and Extended Local Model 
Networks,” Aerospace Science and Technology , Yol. 15, 2011, pp. 33-44. 

19 de Visser, C.C., Mulder, J.A., and Chu, Q.P. “Global Aerodynamic Modeling with Multivariate Splines,” AIAA-2008- 
7500, AIAA Modeling and Simulation Technologies Conference and Exhibit , Honolulu, HI, August 2008. 

20 Barron, A.R., “Predicted Squared Error : A Criterion for Automatic Model Selection,” Self-Organizing Methods in 
Modeling , Farlow, S.J., Ed., Marcel Dekker, Inc., New York, NY, 1984, pp. 87-104. 

21 Morelli, E.A. and DeLoach, R., “Wind Tunnel Database Development using Modem Experiment Design and Multivariate 
Orthogonal Functions,” AIAA Paper 2003-0653, 41 st AIAA Aerospace Sciences Meeting and Exhibit , Reno, NV, January 2003. 

22 Morelli, E.A., “Practical Aspects of the Equation-Error Method for Aircraft Parameter Estimation,” AIAA-2006-6144, 
AIAA Atmospheric Flight Mechanics Conference , Keystone, CO, August 2006. 


Table 1 . Aircraft sensor random error 


root-mean- sq 

uare values 

V, ft/s 

0.05 

a , deg 

0.025 

P , deg 

0.025 

p , deg/s 

0.25 

q , deg/s 

0.25 

r , deg/s 

0.25 
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0.02 

6 , deg 

0.02 

V > deg 

0.02 


0.004 
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0.004 
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0.004 
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